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Abstract 

The evolution of physical and gauge degrees of freedom in the Einstein and 
Yang-Mills theories are separated in a gauge-invariant manner. We show 
that the equations of motion of these theories can always be written in flux- 
conservative first-order symmetric hyperbolic form. This dynamical form 
is ideal for global analysis, analytic approximation methods such as gauge- 
invariant perturbation theory, and numerical solution. 
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One of the prevailing issues facing general relativity, indeed any gauge theory, is the 
separation of physical from gauge degrees of freedom. This conceptual difficulty is encoun- 
tered in generation of solutions of the field equations, proofs of existence and uniqueness of 
solutions, and attempts at quantization. In this letter we present explicitly hyperbolic forms 
]T|] of the Einstein and Yang-Mills equations of motion which clearly display the dynamics of 
these theories without fixing a gauge. (The constraint equations remain elliptic.) The basic 
strategy, which is applicable to any gauge theory, is to take an additional time-derivative of 
the equations of motion, use the constraint equations, and guarantee equivalence to the orig- 
inal theory via appropriate choice of Cauchy data. This completes the program, begun by 
the French, to cast general relativity in 3+1 form by integrating it with the somewhat 
more recent efforts directed towards finding a hyperbolic formulation of general relativity 
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Our hyperbolic formulation preserves complete spatial covariance by means of an arbi- 
trary shift vector. The standard 3+1 treatment 0,|J, is gauge covariant in this sense but 
not hyperbolic. Our formulation does require a condition on the time slicing to deal with 
the time-reparametrization invariance of the theory. 

A hyperbolic formulation of general relativity is valuable for many applications. The 
study of analytic approximations can be given a rigorous foundation. Gauge-invariant per- 



turbation theory |T^] arises naturally as a perturbative reduction of the new equations [|13 



Problems in global analysis, regarding the existence and uniqueness of solutions |T^JT5| , take 



on a new light when viewed with the powerful tool of hyperbolic theory |16| , p!7| ,|8[| . Insights 



into quantum gravity and the problem of time seem likely, given an understanding of the 
precise role of time slicing necessitated by hyperbolicity. 

For numerical relativity, the importance of a hyperbolic formulation cannot be overstated. 
There are many algorithms for solving the hyperbolic equations fluid dynamics which can 
now be applied to general relativity. More fundamentally, the isolation of physical from 
gauge effects means that a numerically generated spacetime can have closer connection with a 
desired astrophysical scenario. In the "Grand Challenge" effort to solve Einstein's equations 



for the inspiral and coalescence of compact binaries — a process expected to be observable 
in gravitational radiation by LIGO and other detectors — one of the difficult problems is 
the treatment of the horizon of the black hole: gauge degrees of freedom can propagate 
faster than light and can thus escape from the black hole. In a hyperbolic formulation 
whose only non-zero speed of propagation is that of light, the horizon again becomes a 
natural physical boundary. The hyperbolic formulation is also ideal for the treatment of 
gravitational radiation in numerically generated spacetimes as it makes manifest a split 
between background and propagating radiation which has long been assumed in approximate 



calculation schemes [18], in the extraction of gravitational radiation waveforms at finite 



radius, and in the imposition of outgoing wave boundary conditions [|T^ . 

We first demonstrate the procedure in the simpler context of Yang-Mills field theory in 
flat spacetime (cf. fl7|l). The Yang-Mills field strength F a ^ u is given in terms of the vector 
gauge potential A a ^ by 

F\ u = d,A\ - d v A\ + f\ c A\A c v , (1) 

where <9 M indicates ordinary partial differentiation in the x M direction in a Minkowski coor- 
dinate frame and f a b c are the structure constants of the Yang-Mills gauge group G. With 
indicating a gauge covariant derivative, the Yang-Mills field equations in the absence of 
sources consist of three equations of motion (for each gauge index value) 

D°F a l0 + D j F a tJ = (2) 

and a constraint 

D j F a 0j = 0. (3) 

The Bianchi identity is 

= D X F\ U + D,F\ X + A,FV (4) 

This is identically satisfied given the definition of F a ^ u and does not need to be separately 
imposed. 



A hyperbolic wave equation for F a io is obtained by taking a covariant time derivative of 
(0) and subtracting a spatial gradient of ([|): 

D D°F a i0 + D D j F a ij - D t D^F a 0j = 0. (5) 

Interchanging the order of covariant differentiations produces (gauge) curvature terms which 
may be combined using the antisymmetry of the structure constants. A covariant divergence 
of the Bianchi identity is then used to give the non-linear wave equation 

D»D,F a l0 + 2f a bc F b ij F c \ = 0. (6) 

The full second-order system of equations consists of the wave equation (§), the constraint 
(0), and the definition of F a ^ in ([!]). [Combining the definition of _F a j in ([!]) with the wave 
equation (||) would produce a third-order hyperbolic equation, with principal part Dd/dt, 
for A a i.] This system is hyperbolic with elliptic conditions for initial data, and its solution 
is unique once Cauchy data have been specified on an initial spacelike hypersurface. The 
Cauchy data consist of an arbitrary gauge potential A a , the pair A a i and F a i0 consistent 
with the constraint (|3p, and D°F a i such that the Yang-Mills equation of motion (0) holds 
on the initial slice. With these data, one can prove the hyperbolic system is equivalent to 
the original Yang-Mills equations, yet no gauge- fixing condition has been imposed |20[| . 

The equations of motion (§) and the definition of F a i0 can be put in flux-conservative 
first-order symmetric hyperbolic form. The magnetic part of the equations, implicit in 
the Bianchi identity, must now be used explicitly. Introducing the derivatives of the field 
strength Gx 1 ^ = D\F a ^ as new variables, one finds 

j^iOri a I r^k^ a r> ra rpb A 77c (n\ 

D (JO + D iJbv — —^J bc-f fx r \v, {< ) 

r->0/^y a r-> siOa ra rpb rpc fo\ 

The unknowns of the first-order hyperbolic system are A a i, F a and Gx""^, and the equa- 
tions consist of the definitions of F a i0 and Gq 11 ^, (|7|) and (||). 

From the first-order form, one sees that A a j, and F a fiu propagate with speed zero: that is, 
they are simply "dragged along" the time axis during the evolution. It is only the derivatives 
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of the field strength that propagate with the speed of light (c = 1). It must be emphasized 
that A a is not a characteristic field and that only the fields which propagate with non-zero 
speed are gauge covariant. A gauge must be chosen to specify A a j, but no gauge- fixing 
condition is required for hyperbolicity. 

A hyperbolic formulation for general relativity can be found by a similar procedure [|TJ. 
(Cf. H where complete spatial gauge covariance is not present because of the choice of a zero 
shift vector.) Consider a globally hyperbolic manifold of topology £ x R with the metric 

ds 2 = -N 2 dt 2 + g^dx 1 + (3 i dt)(dx j + (3 j dt), (9) 

where N is the lapse and f3 l is the shift. Introduce the non-coordinate co- frame, 

6° = dt, 6 i = dx i + ftdt. (10) 

with corresponding dual (convective) derivatives 

d = d/dt - ftd/dx\ di = d/dx\ (11) 

Note that [9 ,9j] = {di(3 k )dk = C 0i k dk, where the C's are the structure functions of the 
co-frame, d9 a = -\C Pl a d^ A 6\ 

The natural time derivative for evolution is 0] 

d = d + f3 k d k -C p = d/dt - Cp, (12) 

where Cp is the Lie derivative in a time slice £ along the shift vector. In combination with 
the lapse as N~ 1 do, this is the derivative with respect to proper time along the normal to 
E, and it always lies inside the light cone, in contrast to d/dt. It has the useful property 
that [do, d/\ =0. The extrinsic curvature Kij of E is given by 

d o9lJ = —2NKij. (13) 

One employs a procedure parallel to that used in Yang-Mills theory. The spatial metric 
gij is analogous to A a i, the shift f3 k to A a , and the extrinsic curvature of E to F a i0 . 
The lapse N is a new feature present in time-reparametrization invariant theories. 
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In four-dimensions, Einstein's theory, R^ v = SnGiT^ — |g M 2,T A A ), leads to six equations 
of motion from Rij, three "momentum constraints" from Roi, and the Hamiltonian constraint 
from G°o = ^(-R°o — R k k)- The hyperbolic form of Einstein's theory is obtained by taking a 
time derivative of the equations of motion and subtracting spatial gradients of the momentum 
constraints, 

doRij — Vj-Roj — Vji?io — &>ij- (14) 

(Barred quantities are defined in the hypersurface £.) 
Expressing (|H]) in a 3+1 decomposition, one finds 

% = NOKij + Jij + Sij, (15) 

where □ = —N~ 1 doN~ 1 do + V fc Vfc is the physical wave operator for arbitrary (3 k . If we 
denote the trace of the extrinsic curvature by H = K k k, then 

J tj = BoiHKij - 2K l k K 3k ) + (N~ 2 d N + H)V l V 1 N 
-2iV- 1 (V fc iV)V (i (iV^ ) ) + 3(V fc iV)V fc ^ (16) 
+N- 1 K ll V k (NV k N) - 2V {l (K 3) k W k N) - NViV 3 H 
+N- l V l V j (N 2 H) - 2NK k {i Rj )k - 2NR kijm K km . 

[where M {ij) = |(M y + M^)} and 

= -N'^iVjidoN + N 2 H). (17) 

For Qij to produce a wave equation, S^ must be equal to a functional involving fewer 
than second derivatives of Kij. This can apparently be accomplished in a number of ways 
and constitutes the imposition of a slicing condition on the spacetime. It is necessary to 
show that the slicing condition can be imposed without spoiling the hyperbolic nature of 
the evolution system. 

A clear and simple slicing condition is the harmonic condition (cf. Q when (3 k = 0) 

d N + N 2 H = 0. (18) 
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(This can easily be generalized by adding an ordinary well behaved function f(t,x) to the 
right hand side.) Imposing fll8| ) for all time amounts to imposing an equation of motion 
for N. The complete system of equations of motion is now the wave equation (|15|) for 



Kij, the harmonic slicing condition (|I8D , and the definition ( |I3| ) of the extrinsic curvature. 
The Cauchy data for the full system, to be given on an initial slice E, are and in- 
consistent with the Hamiltonian and momentum constraints, the lapse N, and d^K^ such 
that the Einstein equations of motion hold on the initial slice. Using the doubly contracted 
Bianchi identity, one can prove |l| that, with these initial data, the hyperbolic system is 
fully equivalent to Einstein's theory. 

Another useful class of slicing conditions arises from choosing if to be a known function 
of spacetime h(t,x). In this case, the lapse function N is determined by solution of the 
time-dependent elliptic problem: 

d h = -V k V k N + N(R + H 2 - g ij Rij). (19) 

In this scheme, the Cauchy data are simply g^ and satisfying the constraints and d^K^ 
from the usual evolution equation. The g^, system is still hyperbolic, but the full set of 
equations is now mixed hyperbolic-elliptic. The shift vector can still be specified arbitrarily. 
Proof of a unique solution proceeds by an iterative method, and equivalence with the usual 
form of Einstein's equations again employs the twice-contracted Bianchi identity [|J . 

The harmonic condition is consistent with the natural slicings of stationary spacetimes. 
For example, suppose one has a spacetime with a timelike Killing vector and a spacelike 
Killing vector proportional to the shift vector: (3 l = In this case the evolution of 

the 3-metric gives dogij = —2^udj^f so H = N^^dif. Eliminating H with the harmonic 
slicing condition yields £ l di(f /N) = 0, which is clearly true for the Kerr geometry in Boyer- 
Lindquist coordinates. It is possible and useful in perturbation theory to have the advantages 
of a specified H and the harmonic slicing condition by choosing the shift vector suitably. 
However, choosing a shift vector in this particular way seems undesirable for numerical 
solution of the full field-equations because spacetimes evolved in this fashion will tend to 



develop coordinate singularities as in the stationary spacetimes mentioned above. 

In the vacuum case, if we introduce, besides N, gij and K^, new variables aj = 
N~ x VjiV — the acceleration of the local Eulerian observers (those at rest in the time slices)— 
its derivatives aoi = N~ 1 dodi and ciji = V well as the derivatives of the extrinsic 

curvature 

doKij = NLij (20) 

and Mkij = V k Kij, one can cast the equations (|13D , (|15|) , flUf ) into complete flux-conservative 
first-order symmetric hyperbolic form 0. The unknowns of the first-order system are g^, 
N, Kij Lij, M ki j, cii, ciji and a 0i , and the equations of the first order system are (|13|), (|T8|), 
(0) and 

a Q Ly - NV k M klJ = NiHLij - Jy), (21) 

<9 il% - NVkLij = N[a k L i:j + 2M k(l m K j)m (22) 
+2i^ m (jMj) fc m — 2K m (iM m j) k 
+2K m (i(K m j)a k + cij)K m k — a m Kj) k )], 

d oai = -N(Ha t + M ik k ), (23) 

d a 3i - NV.aoi = Na k [2M {lj) k - M k ij (24) 
+2a,(iKj) k — a k Kij] + Naja 0i , 

d a 0i - NV k a kl = N[-R k iak + a^H 2 - 2K kl K kl (25) 
+2a k a k + 2a k k ) + 2a k a k l + HM ik k - 2K kl M ikl \, 

where can be found from (|T6|). Notice that the shift is not one of the characteristic fields. 
The form of the first-order system is independent of the choice of /3 fc , though it must be 
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specified for solutions. To complete the reduction to first-order form, the 3-dimensional Rie- 
mann curvature appearing in is expressed in terms of the 3-dimensional Ricci curvature 
using 

Rmijk = 29m[jRk]i + ^9i[kRj]rn + R9rn[k9j]i, (26) 

which in turn is eliminated by substituting 

Rij = Rij + L i:j - HKij + 2K ik K k j + + a^. (27) 

The four-dimensional Ricci curvature is then eliminated using the Einstein equations. 

Note that in spatial dimensions greater than three, the expression for R m ijk involves the 
Weyl tensor, which cannot be eliminated using the Einstein equations. The reduction to 
first order form in these variables is thus blocked. 

One sees that t/y, K^, N, and a* all propagate with zero speed with respect to the 
Eulerian observers: they are dragged along the normal to the foliation by the evolution. Only 
the derivatives of the extrinsic curvature and the derivatives of the acceleration propagate 
with the speed of light. These represent time-dependent tidal forces and can be used to 
form the components of the spacetime Riemann tensor. The only propagating degrees of 
freedom then are curvatures, as one would expect physically. Equivalently, in the second 
order system, the wave equation for can be viewed as defining the notion of radiation as 
distinctly as possible in the context of a nonlinear, curved space field theory. The inherent 
separation of the evolution of the spatial-metric and extrinsic curvature is a natural starting 
point for a formal expansion scheme which is gauge-invariant at each order. 

When linearized around a static background 3-metric, the evolution equations for gij 
and the wave equation for Kij decouple. For instance, with a flat space background, the 
simple wave equation OKy = for the extrinsic curvature is obtained (assuming harmonic 
slicing). The evolution equation for the 3-metric contains no new information about the 
dynamical degrees of freedom. Similarly, one can linearize about static or stationary black 



hole backgrounds, for instance Schwarzschild or Kerr. In the Schwarzschild case, fll5| ) reduces 



directly to scalar wave equations for the even and odd-parity radiation modes. (The 3-metric 
evolution equation is again irrelevant.) Taking an additional time derivative of the scalar 
equations, one recovers the standard results of gauge-invariant perturbation theory [|T2 



The pair of scalar wave equations obtained for each (£, m) multipole combination can be 
matched directly onto a numerically generated interior solution and provides both a gauge- 
invariant radiation extraction method and a clean prescription for outer boundary conditions 
(including backscatter of waves). By refining the assumed exterior background spacetime, 
arbitrary amounts of physical detail can be incorporated by this general method. 

The reasoning we have applied in this paper can be applied to general relativity coupled 
to other fields with well-posed Cauchy problems, as well as to generally covariant and gauge 
theories in the broad sense. One sees that the procedure of taking time derivatives and adding 
further variables can be continued to build a "tower" of equations which, provided suitable 
initial conditions are given, is equivalent to the original theory. Nothing fundamental is 
gained by going beyond the stage at which gauge-invariant equations of motion are obtained, 
as in this paper, but we find the equations that propagate the spacetime Riemann tensor 
components directly aesthetically appealing. By achieving a hyperbolic formulation of a 
gauge theory without gauge-fixing, one has manifestly physical propagation without the 
encumbrance that comes from having to impose particular gauge conditions. The physical 
structure of the theory is revealed with the full gauge symmetry preserved. 

A. A., A. A., and J.W.Y. were supported by National Science Foundation grants PHY- 
9413207 and PHY 93-18152 (DARPA supplemented). 
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